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Dynamical overlap fermions, results with HMC algorithm 
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We present results of a hybrid Monte-Carlo algorithm for dynamical, rhf = 2, four-dimensional QCD with 
overlap fermions. The fermionic force requires careful treatment, when changing topological sectors. The pion 
mass dependence of the topological susceptibility is studied on 6 4 and 12 • 6 3 lattices. The results are transformed 
into physical units. 



1. HMC FOR THE OVERLAP 

It is a great numerical challenge to make dy- 
namical QCD simulations with Ginsparg- Wilson 
fermions, which have exact chiral symmetry at fi- 
nite lattice spacing. We present an exploratory 
study for overlap fermions, which have the Dirac 
operator: 

D = m (l + 7 5 sgn( J ff M /)), 

where Hw is the hermitian Wilson-operator with 
—mo mass pp. The bare mass is introduced by: 

D(m) = to + (1 - to/(2to )) D. 

We use the Zolotarev rational series to approx- 
imate the sign function: 



sgn(x) w e n (x) = x(x 2 + c 2n ) 
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for bi, ci see 0- The evaluation of the inverse of 
the shifted Wilson-matrices is done by multishift 
conjugate gradient algorithm To speed up 
the inversions one can project out the low-lying 
eigenmodes of the Wilson kernel: 

sgn(ff w )«J]sgn(A l )|A l )(A 4 |+Te n (i/ w ), 



where H w \\i) = Aj|Ai) and T = 1 - V\ |A,)(Aj 
projects to the orthogonal subspace. 



We use the HMC algorithm to generate dynam- 
ical configurations for nt = 2 flavors. Thus we 
need the expression of the microcanonical energy: 

n= l -(P 1 P)+S g [U] + ^{D^D)- 1 ^ 



which governs a 
of link variables. 



where ip = (Z? 1 D)~ 
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classical motion in the space 
is the pseudofermion field, 
whereas P is the canonical momenta, which 
is randomized at the beginning of each trajec- 
tory. For the classical motion we need the gauge 
derivative of pseudofermionic action, namely the 
fermion force. 

For the fermion force one needs to invert the 
overlap operator on the 4> field: 

force = — tp^ ■ — t/j, 

dU 

thus one ends up with a nested inversion proce- 
dure, which is the most time consuming part of 
the algorithm. The gauge derivative of the ra- 
tional approximation is straightforwardly calcu- 
lated, furthermore it is easy to include the con- 
tribution of the projected eigenmodes exactly in 
the force 0] . Then one can solve the equations of 
motion with a reversible and area conserving inte- 
gration procedure (usually it is a simple leapfrog). 

2. FERMION FORCE 

With the naive application of the above results 
one faces a further problem. The pseudofermionic 
action has a dicontinuity across the surface, where 
Hw matrix has a zeromode, which means a pres- 
ence of a Dirac-delta type singularity in the force. 
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A finite stepsize integration of the equations of 
motion practically never notice it. This means 
that it drastically violates energy conservation. 
Finally one ends up with a very bad acceptance 
ratio, since the accept/reject step depends ex- 
ponentially on the energy-conservation violation. 
We have observed that already on a 6 4 lattice 
there were no accepted configurations at all. 

These eigenvalue crossings are easily connected 
with the topological sector change, since topolog- 
ical charge Q changes whenever a Xi crosses zero: 

AQ = -i-A[Tr( 75 £>)] = ^A^gii^)] = ±1. 

We cure the problem by considering a classi- 
cal trajectory, when it reaches the zero eigenvalue 
surface. Let the jump in the action on the surface 
be AS, whereas the normalvector of the surface 
is N = (\\dH\y /dU\X). If the momentum is not 
large enough to climb the jump (P 2 < 2 AS), the 
trajectory will be reflected off from the surface, 
its new momentum will be: 

P' =P-2N(N,P). 

If there is enough kinetic energy (P 2 > 2 AS), 
then the system will move into the other topolog- 
ical sector with a new momentum: 



P' =P- N(N, P) + N(N, P) y/l - 2 AS/ (AT, P) 2 . 

One should modify the standard leapfrog al- 
gorithm to accomodate the above momentum 
change. The eigenvalue crossing can happen only 
in the half step, when the gauge fields are up- 
dated. We propose to replace this step with three 
other ones (corrected leapfrog step): 

1. Evolve the gauge fields just to the eigenvalue 
surface (t\ is the time, which is needed for this). 

2. Decide whether reflection or refraction hap- 
pened, and modifiy the momentum according to 
the above rules. 

3. Update the gauge fields in the remaining time 
(t/2 — T\) with the new momenta. 

The algorithm is reversible, preserves area and 
conserves energy only upto 0(t%). It is an im- 
portant goal of future algorithm developments to 
improve on this behaviour. 

The stepsize should be small enough to keep 
track the evolution of eigenvalues, and to avoid of 
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Figure 1. Comparison of the uncorrected (tri- 
angles) and corrected (boxes) leapfrogs. Upper 
part is the energy, lower part is the lowest mode 
of Hw- 



the problem of having 2,3 or more crossing eigen- 
values in one microcanoncial step. 

Fig. compares the evolution of the energy 
and lowest A for the usual and for the corrected 
leapfrog. In the uncorrected case there is a huge 
energy jump at the crossing, whereas in the mod- 
ified case a reflection happens, and TL is much 
better conserved. 

3. NUMERICAL RESULTS 

We have checked our code by a brute force ap- 
proach on 2 4 and 4 4 lattices: we have weighted 
a quenched ensemble with the exactly calculated 
determinant. A complete agreement was found. 
On 4 • 6 3 lattices there is a sharp increase in 
the Polyakov loop at f3 — 5.7, this value of the 
coupling was used for measuring the topology on 
6 4 lattices. The negative quark mass was set to 
too = 1-6, the bare fermion mass was in the range 
to = 0.1. .1.15, the stepsize was r = 0.025 in aver- 
age. At each bare mass roughly 800 trajectories 
were generated. 

The results are plotted on Fig. [21 The up- 
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Figure 2. Topology on N s — 6 lattices. See text. 
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per panel shows the charge history. The average 
topological charge is consistent with zero for the 
total mass range (middle panel). 12 • 6 3 lattices 
were used to fix the scale using ro from Wilson- 
loops. The result is a ~ 0.25 fm for small masses. 
Pion masses were also measured and m\ = Am 
with A w 1 is found in lattice units. Using the 
scale and the pion mass, it is possible to get the 
topological susceptibility in physical units (lower 
panel of Fig. [2J. x( TO ) tends to zero for small 
quark masses. One can compare these results 
with the continuum expectation in the chiral limit 
(solid line of the figure): 

(Q 2 ) />» 



lim x( m ) 

m— >0 



lim 



V 2n f 

Caution is needed, when interpreting the re- 
sults. There are obvious problems such as small 
volume and rough lattice. Furthermore there can 
be doubts about the locality of the theory at these 
rather large couplings [5]. 

A detailed description of the algorithm can be 
found in whereas for an independent study 
see 0. 
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